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Abstract 

y—i We study focusing discrete nonlinear Schrodinger equations and present a new variational 

existence proof for homoclinic standing waves (bright solitons). Our approach relies on the con- 
strained maximization of an energy functional and provides the existence of two one-parameter 
families of waves with unimodal and even profile function for a wide class of nonlinearities. 
Finally, we illustrate our results by numerical simulations. 
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1 Introduction 

The discrete nonlinear Schrodinger equation (DNLS) is one of the most fundamental lattice equa- 
tions and plays a prominent role in the theories of optical waveguides, photorefractive crystals, 
and Bose-Einstein condensates. For an overview on applications and results we refer the reader to 
0^ [KRBOli IEJ031 [Kev09a, IPor09j and the references therein. In this paper we aim in contributing 
to the mathematical theory of DNLS by giving a new variational existence proof for homoclinic 
standing waves. 

(N 
O 

In one space dimension the homogeneous DNLS is given by 

r—i 

> iA j +a(A J+1 + A j ^ 1 -2A J ) + f3A j + *'(\A\ 2 j )A j = 0, (1) 

• i— i 

where j is the discrete space variable, t is the time, and Aj = Aj(t) is the complex- valued amplitude 
or the value of the wave function. The potential function \P is often assumed to be monomial with 
\P'(a;) = x and ^'(x) = x 2 for the cubic and quintic DNLS, respectively. 

For our purposes it is more convenient to simplify the linear terms by means of the gauge 
invariance of DNLS, that is the symmetry under Uj ~» e lip Uj with ip G M. More precisely, with 
Aj ~~> e 1 '^ - 2a '*Aj we readily verify that is equivalent to 



±Aj + a(A j+ i + A,-_i) + *'(|A| 2 )Aj = 0. (2) 
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It is well established that the dynamical properties of ([2]) strongly depend on the sign and strength 
of the coupling parameter a. For convex ^, one usually refers to a > and a < as the focusing 
and defocusing case, respectively. In the anti-continuum limit a — > the DNLS becomes a system 
of uncoupled oscillators, whereas in the continuum limit a —> ±oo the DNLS can be viewed as a 
finite difference approximation of the nonlinear Schrodinger PDE (via the scaling j ~~> y/\l/a\j). 

During the last decades a great deal of attention has been paid to coherent structures such as 
travelling waves and standing waves. Standing waves can be regarded as relative equilibria which 
stem from the gauge invariance. They are special solutions to ^ with Aj(t) = e lcrt Uj(t), where 
the profile u = (uj)- is assumed to take values in R and satisfies 

auj = a(uj + i + Uj-i) + $'(uAuj. (3) 

Standing waves come in different types: Periodic waves (or wave trains) satisfy Uj = iti+jv for 
some periodicity length N < oo. Homoclinic waves (solitons, solitary waves) are localized via 
lim J _ s .-|- 00 Uj = 0, and are hence also breathers. Finally, heteroclinic waves (fronts, kinks) connect 
different asymptotic states. 

The existence of standing wave solutions to ^ has been investigated by several authors using 
rather different methods. Sometimes it is possible to find exact solutions, see jELS85l KRSS05J, 
but in general one needs more sophisticated and robust arguments. A perturbative approach to the 
existence problem was developed by MacKay and Aubry [MA94J IAub97] and relies on continuation 
method. The main idea is to start with a given solution in the anti-continuum limit a = and to 
show that there is a corresponding solution for small a. Continuation arguments have been proven 
powerful for both analytical considerations and numerical simulations and seem to be the preferred 
method in the physics community. 

The main limitation of any continuation methods, however, is the need of an anchor solution 
around which the equation is expanded. As a consequence there is a growing interest in alternative 
existence proofs for standing waves. Example are dynamical systems approaches |QX07[ [ PR05] , or 
variational methods that employ critical point techniques (linking theorems, Nehari manifolds) to 
establish the existence of waves with prescribed frequency a, see |PZ0H IPR 08, ZP09, ZL09J and 
|Pan06[ IPan07l ISZlOj for similar results in DNLS with periodic coefficients. 

1.1 Variational setting 

In this paper we rely on a variational setting, which does not prescribe the frequency but the 
power of a standing wave. More precisely, we obtain homoclinic standing waves as solutions to 
a constrained optimization problem, in which the frequency a is the Lagrangian multiplier. A 
similar idea was used by Weinstein |Wei99j for DNLS with power nonlinearities, but we allow for 
a wider class of nonlinear potentials \£. Moreover, our approach provides more information about 
the shape of standing waves as it guarantees the existence of waves with unimodal and even profile 
u. We also emphasize that, contrary to variational methods with prescribed a, our existence proof 
gives rise to an effective approximation scheme for standing waves. Finally, the restriction to the 
one-dimensional case is not essential but was made for the sake of simplicity. 

In oder to sketch the main idea of our method we introduce an energy functional V and the 
power functional M by 

v(u) = Y,^i u l) + a Yl ("j+i + ^( u ) = Yl u l ( 4 ) 

j 3 j 

Both j\f and V are related to conserved quantities for the Hamiltonian system ([T]). In fact, j\f is 
linked to the gauge invariance by Noether's Theorem, and rearranging the quadratic terms we find 
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V{u) = 2aM(u) -H(u), where 



u(A) = Y,a\A 3+ i-A J \ 2 -n\H 2 ) 

j 

is the Hamiltonian corresponding to ([!]). We readily verify that the standing wave equation ^ is 
equivalent to 

adM(u) = dV(u) (5) 

where d denotes the variational derivative with respect to u. The key observation is that ([5]) can 
be considered as the Euler-Lagrange equation of the optimization problem 

maximize V(u) under the constraint Miu) = g, (6) 

where a plays the role of an Lagrangian multiplier. Notice that ^ is equivalent to minimizing the 
energy T~L subject to prescribed power A/", which is a well established idea in the theory of standing 
waves, see |Wei99j for DNLS and |Pav09l IStu09| for dispersive Hamiltonian PDEs. 

In this paper we refine the optimization problem ^ by considering only those profiles u that 
are non-negative, unimodal, and even. Specifically, we solve the optimization problem 

maximize V(u) under the constraints Af(u) = g and u E C, (7) 

where the convex cone C consists of all profiles u that satisfy u-j = Uj and Uj > %+i > for 
all j > 0. Of course, we then have to show that each solution to the so restricted optimization 
problem satisfies the standing wave equation Q without further multipliers. 

It is known that there exist two possible choices for the index j. In the on-site (or site- centered) 
setting we suppose j E Z, whereas in the inter-site (or bond- centered) setting we choose j E Z+ ~. 
Both settings are equivalent on the level of ^ , and even on the level of ^ with ufl 2 , but lead to 
different results when studying waves with even and unimodal profile u E C. In fact, on-site waves 
attain their maximum in an odd number of points centered around j = (generically only in j = 0), 
whereas the maximum of inter-site waves is realized in an even number of points (generically only 
in j = -\ and j = \). 

1.2 Sketch of the proof and main result 

Due to the lack of strong compactness, is not trivial to show that V attains its maximum on the 
set of interest. A standard strategy would be to employ Lion's concentration compactness principle 
[Lio84] . see also |Wei991 lPav09| . but we argue differently: At first we consider the analogue to ([7]) in 
the space of periodic profiles with Uj = Uj+N- The existence of a maximiser is then granted and the 
invariance properties of the reversed gradient flow for V ensure that the maximizer solves ^ with 
some multiplier a. In particular, there are no multipliers due to the unimodality constraint, and 
so we obtain periodic waves with unimodal and even profile. In the second step we then establish 
the existence of homoclinic waves by passing to the limit N — > oo. To this end we exploit a strict 
maximum condition, which replaces the concentration compactness principle and guarantees that 
the periodic waves are uniformly localized. We mention that approximation by periodic waves 
is also used in jPan06| IPan07] . but in a variational setting that prescribes a and uses critical 
techniques to construct saddle points of an action integral. 

Our main result on standing waves for DNLS can be summarized as follows. 
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Theorem 1. Suppose a > and that \E f satisfies the super-linear growth and regularity conditions 
formulated in Assumption^ Then, in both the on-site and the inter-site setting there exists a one- 
parameter family of homoclinic standing waves (u, a) that is parametrized by g = \\u\\ 2 > g*(a) > 0. 
For each wave, the frequency satisfies a > 2a and the profile u is non-negative, even, unimodal, 
and exponentially decaying. 

We proceed with some remarks concerning the assumptions and assertions of Theorem [T] 

1. The class of admissible includes all convex functions with ^'(O) = 0. In particular, it 
contains power nonlinearities 

= -x 1+r >, V'(x) =x r >, r/> 0, 
V 

but also potentials with saturable derivative such as 

$(x) = x - log(l +x), <&'(x) = 

1 + X 

Notice that the dynamical systems approach by Qin and Xioa |QX07| also allows for arbitrary 
strictly convex ^ and guarantees the existence of a homoclinic wave for all a £ (2a, 2a + h 00 ) 
with hoo = linx^oo *$>'(x). 

2. Theorem [T] also covers some non-convex potentials as for instance 

x 3 r „ , x 2 (3 + x 2 ) 
' l + x 2 ' K J (1 + x 2 ) 2 

3. Our results derived below imply g*(a) — > as a — > 0. Therefore, instead of fixing a and 
choosing g sufficiently large we can alternatively fix g and choose a sufficiently small. We also 
derive explicit criteria for ^, which guarantee that g*(ot) = and reflect well known results 
about the excitation threshold for breathers in ID, see [FKM971 IWeTMl IDZC08] . 

4. Theorem [T] provides the existence of bright solitons for a > 0. Via the staggering transfor- 
mation Uj (—l) J Uj it also implies the existence of standing waves for the defocusing case 
a < 0. These resulting waves have an alternating phase structure and are called bright gap 
solitons, see |Kev09b| . However, the most fundamental standing waves for a < are dark 
solitons which correspond to heteroclinic solutions to ^ and require a different variational 
setting. This is discussed in [HerlOj . 

We do not claim that the constrained maximization of V provides all standing wave solutions to 
Q. It fact, by continuation from the anti-continuum limit it is known that there exist infinitely 
many multipulse solitons, see |Kev09c| . However, the careful spectral analysis by Pelinovsky et 
al. [PKF05J indicates that most of these multipulse solitons are unstable. This is true even for 
unimodal inter-site waves, but at least unimodal on-site waves can be expected to be stable. We 
also emphasize that maximisers of the optimization problem ^ can be shown to be orbitally 
stable Wei86, Wei99j, for instance using the Grillakis-Shatah-Strauss theory, see |Pav09] and the 
references therein. Due to the additional shape constraint u G C, however, this stability result 
does not apply directly to the waves provided by Theorem [T] because it is not clear that the 
optimization problems ([6| and ([7]) have the same solution. We conjecture that this is indeed the 
case but a rigorous proof is still missing. 

This paper is organized as follows. In ^2] we formulate our assumptions on ^ and introduce some 
notations. We also introduce the strict maximum condition and investigate the maximum of V on 
bounded subsets of i 2 . In ^3] we then proof the existence of periodic waves and pass to the limit 
N — > oo. Finally, we present some numerical simulations in £j4j 
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2 Setting and properties of the energy landscape 

In this paper we always assume a > and rely on the following standing assumptions on 

Assumption 2. The potential Vl/ is continuously differentiable on [0, oo) and has the following 
properties. 

1. ^ is normalized by ^'(0) = \j/(0) = 0. 

2. VP is grows super-linearly, that means x^'(x) > *$>(x) > for all x > 0. 

3. is non-degenerate, i.e., ^>(x) > for x > 0. 
Remark 3. 

1. Assumption [2] implies 

< x~ 1 ^(x) < &(x) ^> 0, *(As) > A*(x) Vx>0, A>1, <S(x) oo. (8) 

2. Each convex and normalized potential \£ grows super-linearly since the function 0(a?) = 
x^'(x) - $f(x) is non-negative due to 9(0) = and 9'(x) = x^"(x) > 0. 

3. Assumption [2] is invariant under scalings 

$( x ) := a^(6x 1+c ), 

where a,b,c> are given constants. This means ^/ satisfies Assumption [2] if and only if ^ 
does so. 

2.1 Notations 

In what follows we consider real-valued sequences u = (uj)j eJ , called profiles, where the index set 
is given by J = Z and J = Z + | in the on-site and inter-site setting, respectively. Under the 
periodicity condition Uj = uj + n for all j and some N, each profile u is uniquely determined by its 
values on a periodicity cell Zjy. We choose 

Z 2A / = {—M + 1, .., -1, 0, 1, .., M}, Z 2M+1 = {—M, -1, 0, 1, .., M}, 

in the on-site setting, and 

Z 2M = {— M + 1, -1 i, M - §}, Z 2M +i = {—M + 1, -i i, M + §} 

in the inter-site setting. The symmetrized periodicity cell is abbreviated by Zn = Z^ n (— Zn) 
and we readily verify, see Figure [1} that 

(iV-l)/2<maxZjv < N/2, #{i G Z N : \i\ < \j\ } = 2 \j\ + 1 VjGZjv- (9) 

For our purposes it is convenient to identify the non-periodic case with N = 00, so we write 
Zqq = Z&Q = J. In view of Q and (§ we introduce the discrete function spaces 

= ju = (uj)j eJ : % = uj + n for all j G J, ^ it| < 00 1 = R N , 

jez N 

£l = {u= ( Uj ) jeJ : Yl u ) < °°} = ^( Z )' 
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Figure 1: Examples for periodic profiles that are non-negative, even, and unimodal. The values in 
the periodicity cell are plotted in Black. 

for N < oo and N = oo, respectively. For both finite and infinite N we denote by 

INI 2 = Yl n i' v ) = Yl u i v i> 

jeZjv j&z N 

the norm and scalar product in t 2 N , so the sphere of radius ^/~g reads 



'Sn, q = {u g 

We also define the energy functional "Pat by 



f2 



g}. 



V N {u) = aC N (u) + W N (u), C N {u) = Uj(u j+ i + Uj_i), 

j&z N 

and denote by 



jez N 



Cn = {u : Uj = u-j > V j G Zjy} n {u : Uj-± > Uj V 1 < j G Z/y} 

the set of all profiles u that are non- negative, even and unimodal on the periodicity cell Zjy. 
From these definitions we readily draw the following conclusions. 

Lemma 4. For both finite and infinite N the following assertions are satisfied. 

1. Vn is Gateaux- differentiable on l 2 N with derivative &P(u)j = 2a(uj + \ + ity-i) + 2^'(u 2 )uj. 

2. Cn is a positive convex cone and closed under weak convergence in £ 2 N . 



3. By Q we have Uj < ||it||/y2 \j\ + 1 for all u G Cn and j G Zjy. 

The core of our variational existence proof for standing waves is the optimization problem 

maximize Vn on the set Cn n Sn, e - (10) 

For finite N, the existence of maximizers follows from simple compactness arguments and the 
invariance properties of the gradient flow of Vn imply that each maximizer is in fact a periodic 
standing wave. For infinite N, however, we lack compactness and construct homoclinic maximizers 
of V oo as limits of periodic maximizers of Vn- 
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2.2 Energy maxima on bounded sets 

In this section we investigate the function 

1 

T N (a, q) = — sup{V N (u) : u G C N n S Nj g } , 
ag 

and show that the strict maximum condition 

TooCa, f?) > 2 (11) 
holds if one of the following conditions is satisfied: 
(^41) a is sufficiently small, 
(A2) g is sufficiently large, 

(A3) ^(x) ~ cx 1+r ' for all < x < 1 and < i] < 2 and some c > 0. 



The strict maximum condition (11) appears naturally in our existence proof for homoclinic waves, 
see £3.2, and guarantees that the influence of the nonlinearity VP is strong enough in the limit 
N — > co. Specifically, if (11) is satisfied, then the optimization problem (10) has a solution for 
N = co, and each minimizer is a homoclinic wave with unimodal profile. Moreover, if (11) is 
violated, then we have T OQ (a, g) = 2 = sup Mgi5oo £oo(u), and this implies that (10) has no solution 
for N = co. 

We also emphasize that (11) is equivalent to 



inf jft(u) ; u G £ 2 (Z), Af(u) = g} < 0, 



which is precisely the condition used in [Wei99j to prove the existence of homoclinic standing waves 
for DNLS with power nonlinearity = cx 1+n . Moreover, condition (A3) implies that there is 
no excitation threshold for power nonlinearities with < rj < 2, which means there exist waves 
with arbitrary small g. This result is again in line with the findings from |Wei99j. 

We now summarize some elementary properties of the function Tjy. 

Lemma 5. Tjy is decreasing in a and increasing in g for both finite and infinite N. 

Proof. The monotonicity with respect to a is obvious. Towards the monotonicity in g let u G i\ 
and A > 1 be fixed. We have Cn(Xu) = X 2 Cn(u), and 

-f x W N (Xu) = £ 2*'(A 2 «f)Au 2 > £ -^(X 2 u 2 ) = - x W N (Xu) 

jez N j&z N 

implies Wn(Xu) > A 2 WV(u). Consequently, for ||u[| ^ and A > we find 

T N (a, X 2 \\u\\ 2 ) > -^^(C N (Xu) + ^W N (Xu)) > J—(£ N (.u) + JW«))> 

and taking the supremum over u gives T^{a, A 2 ||ii|| 2 ) > TpfCu, ||w|| 2 )- Q 

Remark 6. For power nonlinearities ty(x) = cx 1+v with c, r] > we have T/v(a, A 2 £>) = 
jv(A~ 2r? a, and hence 

T N (a, g) = T N (ag~\ l) = T N (l, a'^g) 
for all a, g, A > and both N < oo and iV = co. 
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For finite N, the constant profile with uj = \J g/N for all j belongs to Cn DSn,q, and using ^ 
we readily verify that 



V n {^/q~Jn) =2ag + g—^(g/N) 2ag. (12) 



In particular, due to *$>(g/N) > we have 



T N (a, g) >V n {Vq/N) >2 (13) 
for all N < oo, a > 0, g > 0. The case N = co is a bit more delicate. 

Lemma 7. We have Too (a, £>) > 2 wt/i strict inequality provided that (Al) or (A2) is satisfied. 
Proof. We start with the proof in the on-site setting and define u m 6 n 5oo,£ by 

_ \/g f 1 for |j| < m, 
Mm,J ~ ^2^ + 1 [0 for |j| > m. 

A direct calculation shows 

2m ( g \ 

Coo(u m ) = - —g, Voo(u m ) = (2m + 1)V — , 

2m + 1 \2m + 1/ 

and hence 

aCoo{u m ) + Voo{u m ) 



Too(a, g) > 



ag 

1 g + 12m+l^_^_\ >2 . 1 



£> 2m + 1 a £ \ 2m + 1 / £> 2m + 1 



(14) 



This gives Too (a, g) > 2 by passing to the limit m — > co and using 

We have now two possibilities to show T OQ (a, g) > 2. First, for given g and m we can choose a 
small. Second, for fixed a and any g > 1 we choose m = m(g) such that 2m — 1 < g < 2m + 1, so 
(14) gives 

Too(q, g)>2-- + - 
g a 

for some c > 0, and the claim follows with g — > co. 
Finally, in the inter-site setting we define u m by 



y/Q j 1 for \j\ <m-\ 



V 2m I for \j\ > m 



and argue analogously. □ 

We now show that the strict maximum conditions is always satisfied if \£ is sufficiently nonlinear 
for x ~ 0. The key idea in our proof is also used in |Wei99] to disprove the existence of an excitation 
thresholds for power nonlinearities with < rj < 2. 

Lemma 8. (A3) implies T OQ (a, g) > 2. 
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Proof. We present the arguments for the on-site setting; the proof for the inter-site setting is 
similar. For fixed a > 0, g > but arbitrary < ( <C 1 we consider u^ G Coo n Soo, e defined by 



u (jj = v / £»tanhCexp(-C|i|), 
which has small amplitudes due to tanh£ = £ + 0(C 3 )- By direct computations we find 

£ooK) = ^ = ,(2-c 2 + o(C 4 )), 

as well as 

WU«c) = (l + °(l)) £ cn/^ = (1 + o(l)) c( g tanh C ) 1+ ^ + ^ \ - ^ ± ^ 

1 - exp ( - 2C(1 + ??)) 

= (1 + «(1))^ 1+ ^, 
where o(l) means arbitrary small for small £. We now conclude that 

Too(a, p) > ^K) = 2 + (1 + o(l)) (^^yC " C 2 

and the claim follows from choosing £ sufficiently small. □ 

We finally establish the convergence T/v(a, g) — > T OQ (a, g) as N — > oo. To this end we denote 
elements of l N by Un = («j\r, j)j and introduce two operators Rn and i£jv which act on profiles u 
as follows. Rnu is the restriction of u to the symmetrized periodicity cell Zn, i.e., 

(Rnu) ={ Uj ioT ^ e ^ N ' H5) 
J I otherwise, 

and i?7v^ is defined as the periodic continuation of Rnu. Notice that we allow for small embedding 
errors. In fact, for non-symmetric periodicity cells with Zn ^ Zn we have Enun ^ un for 
un G £ N , but the embedding error is small due to the decay estimate from Lemma |4j 
We readily verify that En and Rn are linear operators with 

Rn ■ Cn — >■ Coo, En ■ Coo — >■ Cn- 

and ||i?jv"itjv|| < H^aHI an d H-^JV^ooll < ll u oo||- 
Lemma 9. We have 

\V N (U N ) ~ Voo{R N U N )\ < OQ\u N \\/VN), Vn(EnUoo) 7>oo(«oo) 

for all un G Cn and Uoo G Coo- 
Proof. The first claim is a consequence of ^ and Lemma |4j The second one follows from RnUoo — > 
Uoo as N — > oo. □ 

We are now able to prove the convergence of maxima. 
Lemma 10. We have Tn(o, g) — > Too(a, g) as N — > oo. 
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Proof. For fixed Uoo G Coo <~^<Soo,g, Lemma [5] and Lemma [9] provide 

Voo(uoo) < Vn(EnUod) + o(l) < agT N (a, ||^ArUoo|| 2 ) + o(l) 
< agT N (a, g) + o(l), 

with o(l) being arbitrary small for large N. Passing to the limit N — > 00 and taking the supremum 
over Uoo gives Too(a, g) < timinfjv_> DO T^{a, g). Towards the reverse estimate, we denote the 
maximizer of Vn in C-N H <Sjv, g by un, and thanks to Lemma [9] we find 

agT N (a, g) = V N {u N ) < Voo(E N u N ) + o(l) < a^r oo (a, g) + o(l), 

which implies limsup^v^^ T/v(a, g) < Too(a, g). □ 



3 Existence of standing waves 
3.1 Periodic waves 

We now prove the existence of periodic waves with ./V < 00. To this end we consider the reversed 
gradient flow for Vn under the constraint ||u|| 2 = g, that is 

£u = F N (u), F N (u) = dV N (u)-a N (u)u, M^ ^y]'^ , (16) 

where r > denotes the flow time, and cat is a (dynamical) Lagrangian multiplier. Constrained 
gradient flows are well known in the context of ground states for Schrodinger equations and some- 
times referred to as imaginary time methods, see [B D04| IYL07| and references therein. 



Lemma 11. The gradient flow (16) has the following properties : 
1. Sn,q is an invariant set. 



2. Vn(u) is strictly increasing on each non-stationary trajectory. 

3. Each stationary point u £ i 2 N of (16) is a standing wave with frequency a = o~jy(u). 

dr I 



Proof. The definition of a implies -^-||ii|| 2 = (-Fjv(it), u) = 0, so ||it[| is conserved under the flow. 



By a direct computation we find 

£V N (u) = {dV N (u), F N (u)} = J_(||^( u )|| 2 ||u||2 _ \(&p N (u), n)| 2 ) > 

with -^:Vn(u) = if and only if u and &Pj^{v) are collinear, i.e., 8Vn{u) = ctn(u)u. □ 
The next result provides a key ingredient for our existence result. 



Lemma 12. Cn is invariant under the gradient flow (16). 

Proof. Within this proof let u(0) G Cn be some fixed initial data and denote by u(r) G i 2 N with 
r > the corresponding forward in time solution to (16). We first note that 

£n = {u : Uj = U-j V j G Zjy}, 

the set of all even profiles, is a closed linear subspace of £ 2 N and moreover invariant under Fn- We 
thus conclude that u(0) G £n implies u(r) G £n for all r. In order to show the remaining assertions 
we proceed as usual and consider the perturbed initial value problem 



d- 



v = F N (v)+e, v(0) = u(0) + e (17) 
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where the perturbation e = is supposed to be an inner point of Cat with respect to the induced 
topology of £n- First suppose for contradiction that there exist To > and jo £ Zn such that 

v jo( T o) = and Vj(r) > for all j £ Zn and < r < tq. 

This implies -§pVj (jo) < 0, so (17) provides 

-£fV jo (T ) = F NJo (v(t )) + e jo = a(v jo+1 (T ) + u io _i(r )) + e jo > e jo > 0, 

which is the desired contradiction. Secondly, assuming 



i(tq) = Vj (ro) and Vj-\(r) > Vj for all 1 < j < maxZ^ and < r < tq, 



we find a contradiction by 

> ^ru io -i(r ) - -£:V jo (T ) = a(v jo - 2 + v jo - v jo+1 - u io _i) + e io _i - e io > e jo _i - e io > 0. 

In conclusion, we have shown v (r) E Cat for all r > and all perturbations e, and n(r) £ Cat 
follows by ||e|| — » 0. □ 

Theorem 13. For each g > there exists a standing wave (u, a) with u S Cm n 5 at, ^ and ag > 
Vn(u) = agT^^a, g) > 2ag. 



Proof. Due to N < oo the set Sn, g H Cat is compact in 



f^, and hence there exist a maximizer u 



for "Pat on this set. Lemma 11 combined with Lemma 12 then implies that u is stationary point 
of (16) and solves the standing wave equations Q with frequency a = o~n(u). Finally, testing ^ 
with u gives 



a N (u)g = aC N {u) + ^ ^ > a£jv(«) + ^ = ^jv( 



u 



due to Assumption [2j and ( 13 ) completes the proof. 
We conclude this section with two remarks. 



□ 



1. Theorem 13 provides the existence of two families of periodic waves as it holds in both the 



on-site and the inter-site setting. 



2. Theorem 
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does not exclude that the maximizer is equal to the constant profile \J g/N. 
However, Lemma [7J combined with Lemma |10| and (12), ensures that the profile is non- 
constant for large N, provided that (Al), (A2), or (^43) is satisfied. 

3.2 Homoclinic waves 



In this section we prove that the period waves from Theorem 13 converge to homoclinic waves 



provided that the strict maximum condition (11) is satisfied. To this end we fix g > and consider 
a sequence of profiles (uat)at c such that 



1. un is the image of a maximizer of Vn ° n Cn H 5at, g under the restriction map Rn from ( 15 ), 

2. cjat is the corresponding frequency. 
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According to these definitions and Lemma [4] we have 

< u NJ = y/g/(2 \j\ + l) V j, u NJ = V \j\ >(N + l)/2. (18) 

Moreover, Lemma [TU] and Theorem |13| provide 

vnunj = a(u N>j+ i + unj-i) + ®'(u Njj )u N j for all j, N with \j\ < (N - l)/2, (19) 
as well as 

°NQ > VooiuN) + o(l) = agT^g, a) + o(l). (20) 
We next show by using the strict maximum condition that the profiles un are localized. 



Lemma 14. Suppose that (11) is satisfied. Then, 

a = lim inf o~n > 2a 

iV-»oc 

and there exist two positive constants C and d such that 

unj < Cexp (-d \j\) 

holds for all j, N. 



(21) 



Proof. The first claim follows from ( 20 ) . Now choose and j± > 1 such that 



2a < a* < a, 



sup 

0<x<£?/(2>+1) 



Combining this with ([18]), (19), and unj > unj+i gives 
(cr* - a)u N)j < (a N - ^'(u 2 Nj ) - a)u N) 



< 



O-NUNJ - ^'( u N,j) u N,j - CtUN,j+l < aUN, 



and hence UN.j < 3 *-\Tq with k = - a _ a < 1 for all j with < j < N/2 — 2. Finally, (21 ) follows 
with d = — In k and C sufficiently large. □ 



Corollary 15. Suppose that a and g are chosen such that (11) is satisfied. Then, the sequence 
( u n)n * s strongly compact and each accumulation point u^ satisfies Uqo G C, ||iioo|| 2 = Q, and 
"Poo(^oo) = ctgT OQ (a, g). Moreover, u^ decays exponentially and is a standing wave with frequency 
(Too > 2a. 

Proof. By compactness we can extract a (not relabelled) subsequence such that ujy — 1 u^ G C in 



, and this yields the pointwise convergence un 



a. 



oo, J 



for all j. The uniform tail estimate 



(21) then implies \\u c 



g and that Uqo decays exponentially for j — > ±oo. We conclude that 



u 7v — ^ Uoo strongly in l 2 ^ as well as Vcx>(uoo) = ^^■n~>oo'P 00 (un) = agT(a, g) > 2ag, where we 
used (20). Moreover, (19) with fixed j and (J20j imply cj/v — > for some > 2a, and exploiting 
(19) for all j we infer that is a standing wave with frequency a^. □ 

We have now finished the existence proof for standing waves. In particular, Theorem [TJ follows 



from Lemma [7] Lemma [8] Theorem 13] and Corollary 15 We finally recall that for given a > 
the condition (A3) on ^ implies (11) for all g > 0, and hence the existence of homoclinic waves 
with arbitrary small energy. 
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4 Numerical examples 

In this section we illustrate our analytical results by numerical simulations of standing waves with 
N < oo. To this end we define a map / : <Sjv, g —> <Sn, g by 



I(u) = ^fg- 



u + tFn{u) 



\u + tF n {u)\Y 

where r > is sufficiently small, and construct standing waves as limits of the iteration 

uo = itini, u k+ i = I(u k ). (22) 

This scheme preserves the constraint ||u|| 2 = g exactly, and is a discrete analogue to the gradient 
flow (16) due to I(u) = u + tF^(u) + o(r 2 ). In order to compute a good guess for the initial profile 



Mini we start with the ansatz 



ki + n 2Xj + k 3 (1 + cos (irj/N)) + k 4 exp ( - 20(j/N) 2 ) , X j 



1 for |j| < 1, 
otherwise, 



where the parameters k% are positive and coupled by the constraint ||ttim|| 2 = Q- Then we sample 
the set of all admissible parameters by about 100 points, and solve the discrete maximization 
Vn{u\n\{Ki)) — > max to find the optimal values for the parameters Kj. 



Initial profile 



Profile after 30 steps 



Residual after 30 steps 



1-1.5 F 



+ 1.: 



ff).0 l»»* '«»» 

-12 12 

Initial profile 

fl.5F " I 



+6. 10" 



+n.n ^»»*»»»* »»»»»»»«! -6.m- 5 l *»»*. . .*««* 



-12 12 

Profile after 30 steps 



-12 12 

Residual after 30 steps 



+ 1.7 



+2. icr 



-12 



12 



+0.0 



-12 



MsJ -2.10- 5 ^ 
12 -12 
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Figure 2: Periodic on-site (top row, N = 25) and inter-site (bottom row, N = 24) waves for the 



data from (23) with r = 1. 



Profile for g=2.38 

+0.24F jf-K. 



+0.24B?^ 

-20 



20 



Profile for g=2A 



P N -2ag versus g 



hi. 50 



fO.OO 



+6. 00F 



-20 



+0.00h 

20 2 2.4 



Figure 3: Periodic on-site waves for the data from (24) with r = 1 and iV = 41. 



Figure [2] shows numerical results for 

ty(x) = x — arctanx, a = 1, g = 10, 



(23) 
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and provides evidence that the algorithm (22) produces a standing wave in both the on-site and 
inter-site setting. Notice that \P satisfies Assumption [2] and is saturable due to lirnj^oo ty'(x) = 1. 
A second example concerns 



exp [x ) 



2 X 



1, 



a 



1, g€[2,3]. 



(24) 



and is shown in Figure [3| The simulations indicate that the periodic on-site waves for ( 24 ) exhibit 
quite different properties for small and large values of g: For g < 2.38 we observe that almost all 
lattice sites are excited and that the profile has small amplitude and is almost constant. Moreover, 
the energy Pjy = ^n^un) is only slightly larger than 2ag. For g > 2.4, however, the profile is 
strongly localized and Pjy is considerably larger than 2ag. We therefore expect that for N — > 00 
the periodic waves converge pointwise to zero and a non-trivial homoclinic wave for small and 
large g, respectively. Notice that this is in accordance with our theoretical results: Since we have 
^(x) ~ for small x, the existence of homoclinic waves is guaranteed only for sufficiently large 
g. 



Profile for N=10, a= 0.5 



+0.99F 



+0.0( U ,»,«,« »,»,» 

-5 5 

Profile for N=10, a= 2.0 

+0.45F »~ • 



+0.45Lt 

-5 



Profile for N=50, a= 0.5 



fO.99 



P N -2ag versus N 



fO.OO 
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-25 25 

Profile for N=50, a= 2.0 



10 30 50 

P N -2ag versus N 



f0.20 



f0.20 



f0.02R 



4-0.01" 



fO.OOL 



-25 



25 



10 



30 



50 



Figure 4: Periodic inter-site waves for the data from (25) with r = 1 and several values for N. Top 



and bottom row correspond to a = 0.5 and a = 2.0, respectively. 

A similar phenomenon can be observed in Figure HI which illustrates the limit TV — > 00 for 



a G {1/2, 2}, g = 2. 



(25) 



For sufficiently large a we have P/v - > 2ag as N — > 00 and the periodic waves converge (weakly in 
I 2 ) to zero. If a is sufficiently small, however, we have lini/v-!.oo fjv > 2ag and the periodic waves 
converge (strongly in £ 2 ) to a non-trivial homoclinic wave. 
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